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We report an observation of a commensurate-incommensurate phase transition in a Dzyaloshinskii- 
Moriya spiral magnet Ba2CuGe2C>7. The transition is induced by applying a magnetic field in 
the plane of spin rotation. In this experiment we have direct control over the strength of the 
commensurate potential, while the preferred incommensurate period of the spin system remains 
unchanged. Experimental results for the period of the soliton lattice and bulk magnetization as a 
function of external magnetic field are in quantitative agreement with theory. 



Studies of Commensurate-Incommensurate (CI) phase 
transitions have a long history, dating back to pioneering 
works of Frenkel and Kontorova jj] and Frank and van 
der Merwe. 0. Since then CI transitions were discovered 
and studied in a number of such seemingly unrelated sys- 
tems as noble gas monolayers adsorbed on graphite sur- 
face HQ, charge density wave materials H|), an d rare- 
earth magnets [Q (For comprehensive reviews see for ex- 
ample Ref. Q). As a rule, CI transitions result from 
a competition between two distinct terms in the Hamil- 
tonian that have different "built-in" spatial periodicities 
and are often referred to as potential and elastic energy, 
respectively. The potential energy by definition favors 
a structure commensurate with the crystal lattice. The 
elastic term is intrinsic to the system where the transi- 
tion occurs, and has a different "natural" built-in period. 
In many known realizations of CI, such as adsorbed gas 
monolayers, it is the period set by the elastic term that 
can be varied in an experiment to drive the transition, 
whereas both the strength and the period of the poten- 
tial remain constant. In other systems, such as rare- 
earth magnets, both the elastic term (exchange coupling 
between spins) and the potential (magnetic anisotropy) 
can be changed, but only indirectly, by varying the tem- 
perature. 

From the very start it was clear that in its purest form a 
CI transition may be driven by the change of the strength 
of the potential alone, with the two "built-in" periods 
remaining constant. An elegant experimental realization 
of this type of CI was first proposed by Dzyaloshinskii 
||, who considered an incommensurate spiral magnetic 
structure in a magnetic field applied in the plane of rota- 
tion of spins. In this model incommensurability is intrin- 
sic to the spin system and results from spin-spin inter- 
actions. The role of the potential is played solely by an 



external magnetic field H, that favors a commensurate 
spin-flop state. In a real magnetic material with this 
type of CI transition an experimentalist would have a 
convenient handle on the strength of the potential term, 
adjusting it by simply changing the magnetic field. Such 
systems are not easy to find. In most spiral magnets, 
cubic MnSi (k]] 
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e.g., cubic MnSi |iL]J and FeGe 
sufficient to realign the spin plane perpendicular to the 
field direction. While RbMnBr3 [jl2j may actually be 
one material where this does not occur, to our knowl- 
edge, the most crucial quantities, namely the incommen- 
surability parameter £ and magnetization M, have not 
been measured as a function of the external field in any 
"clean" realization of Dzyaloshinskii's model to date. In 
the present paper we report the first direct experimen- 
tal observation of a Dzyaloshinskii-type field-induced CI 
transition in Ba2CuGe207, also presenting experimental 
data for ((H) and M(H). On the theoretical side we 
go beyond a qualitative analysis of the critical properties 
(close to the phase transition) , as was previously done by 
Dzyaloshinskii. Dealing with this particular system, we 
construct an exactly solvable model and derive exact re- 
sults that are in quantitative agreement with experiment 
throughout the entire phase diagram. 
The 

structural and magnetic properties of Ba2CuGe207 are 
discussed in detail in Ref. |L3) . In the layered tetragonal 
crystal structure the magnetic Cu 2+ sites form a square 
lattice with nearest-neighbor (nn) distances of 6 A along 
the (1,1,0) and (1,-1,0) directions, respectively. The low- 
temperature (Tjy = 3.26 K) magnetic phase is a weak dis- 
tortion of the Nccl spin arrangement, with all spins con- 
fined to the (1,-1,0) plane and the staggered magneti- 
zation slowly rotating upon translation along the (1, 1,0) 
direction. The propagation vector is (1 + £, C, 0), where 
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C = 0.027. Only the nearest-neighbor in-plane antifer- 
romagnetic exchange constant is significant and is equal 
to J a b — 0.48 meV. The coupling between adjacent Cu- 
planes is ferromagnetic, with | J c |/| J ao \ « 1/37. We have 
previously suggestedthat the incommensurate structure 
is a result of Dzyaloshinskii-Moriya (DM) an- 
tisymmetric exchange interactions. The corresponding 
term in the Hamiltonian can be written as D(Si x S2), 
where x denotes the vector product and D is a vector 
associated with the oriented bond between the two inter- 
acting spins. 

From the symmetry properties of the lattice we de- 
duce that for a bond between two nn Cu sites along 
the (1,1,0) direction, the only allowed components of 
D are those along (1,-1,0) and (0,0,1), respectively 
(Fig. [I]). The (1,-1, 0)-component does not change sign 
from one bond to the next, while the (0, 0, 1) compo- 
nent is sign-alternating. It is the uniform component 
that is liable for the incommensurate distortion of the 
Neel structure, and for the rest of the paper we shall ig- 
nore the oscillating component, assuming that D is in the 
(a, b) plane. The interaction energy is minimized when 
all spins are perpendicular to D. The total exchange 
energy of the pair of nearest-neighbor spins is given by 
2 J ao cos (f> + Dsm(f> — y/ 4 J 2 b + D 2 cos(0 — a), where (j> 
is the angle between spins, and a = arctanZ3/2J Q ;,. The 
energy is a minimum at <fi = tt + a. The classical ground 
state is therefore a spin-spiral, with all spins in the (1,- 
1,0) plane, and the angle between subsequent spins equal 
to 7r + a. The angle a is related to the propagation vec- 
tor £ by a — 2tt( rj 10°. Precisely this spin configuration 
was found in the initial zero-field neutron diffraction ex- 
periments p3[ . 

The central experimental result of this paper is the 
observation of a CI transition in Ba2CuGe207, induced 
by a magnetic field applied along the (0,0,1) direc- 
tion. Single-crystal magnetization measurements were 
performed using a conventional DC-squid magnetome- 
ter in the temperature range 2-300 K. Experimental 
= jM for T = 2 K are shown in Fig. |. For 
H||a no anomalies are observed. In contrast, when the 
field is applied along the c axis, a distinct feature is seen 
around H = 2 T and indicates the presence of a mag- 
netic phase transition. Neutron diffraction experiments 
were carried out on the H9 (cold beam) and H4M (ther- 
mal beam) 3-axis spectrometers at the High Flux Beam 
Reactor (HFBR) at Brookhaven National Laboratory on 
a « 4 x 4 x 4 mm 3 single-crystal sample, in the tem- 
perature range 1.3-5 K and magnetic fields up to 6.5 T, 
applied along the c-axis of the crystal [1(|. In zero field, 
at T = 2.4 K< T N , elastic scans along the (1 + (,(,Q) 
direction show magnetic Bragg reflections centered at an 
incommensurate position ( = ±0.0273 (Ref. jl3), Fig. 
3b). As the magnetic field increases, the peak moves 
in closer to the antiferromagnetic zone-center at (1,0,0) 
[Fig. |, insert (a)]. At H > H c « 2.3 T the satellites at 




(1±C, ±£, 0) are no longer observed, but are replaced by a 
single peak at the C-point (1, 0, 0) [Fig. ||, insert (b)]. The 
magnetic structure thus becomes commensurate. The ex- 
perimental field-dependence of the propagation vector £ 
is shown in the main panel of Fig. |^. 

To quantitatively describe the field-dependent be- 
haviour, we follow the approach of Dzyaloshinskii [[)). We 
assume that the vector of local staggered magnetization 
at point r remains in the (1,-1,0) plane and forms an- 
gle 8(r) with respect to the c-axis. The free energy per 
Cu-plane in the continuous limit is then given by: 



F = 



Here the axes x, y and z run along the (1, 1,0), (1, —1,0) 
and (0,0,1) directions, respectively, and r = (x,y,z). 
The first term is the total elastic energy of exchange in- 
teractions, and favors a spiral structure of period 2zeA. In 
Eq. (Q) A is the in-plane nn Cu-Cu distance, and p s is the 
in-plane spin stiffness, that for classical spins at zero tem- 
perature is given by p s {0) = \J J 2 b + D 2 S 2 . 7 is the spin 
stiffness anisotropy defined by 7(A C /A) 2 = | J c |/| J a b\ ~ 
1/37 for Ba2CuGe207. The second term represents the 
Zeeman energy. X\\{T) and X-i-i^) are defined as mag- 
netic susceptibilities with respect to fields that rotate 
along with the spiral structure and are parallel or per- 
pendicular to the local staggered magnetization, respec- 
tively. In the paramegnetic phase X\\(T) — x_l(T), while 
at T — the classical result is x^l(O) = (gps) 2 /(8J Q hA 2 ), 
and xy(0) = 0. The equilibrium spin configuration 
should minimize the free energy (pi), and therefore satisfy 
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where T = (p s /H 2 ( X ± ~ X\\)) 1/2 - 

Expression (J2j> has the form of the sine-Gordon equa- 
tion, that is central to describing CI transitions in many 
systems, and its "soliton lattice" solutions are well- 
known: 



0(x) = am(x//3T,/3), 



(3) 



where am{x, (3) is the Jacobi elliptic function of modulus 
j3. Analogues of Eqs. (|l|) and (^) were derived in Ref. [||. 
To obtain exact results for ((H) and M(H), however, we 
make one additional crucial step. For each value of H 
of all valid solutions, labeled by f3, one has to choose the 
one that indeed corresponds to the global minimum of the 
free energy. This is done by substituting Eq. (^) into Eq. 
(0) and minimizing the resulting expression with respect 
to p. After some algebra we find that F is minimized 
when 
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where E(f3) is the elliptic integral of the second kind. H c 
is the critical field at which the CI transition occurs |j] . 
Indeed, for H > H c the spin structure is given by the 
soliton-free solution 8(r) = ir/2, which corresponds to a 
commensurate spin-flop phase, as visualized in Fig. ^(a). 
The staggered magnetization in this case is parallel to 
the x-axis and the spins are slightly tilted in the direc- 
tion of the field. In the limit H = one has (3 — > 
and 6(x) = |-f, which corresponds to an unperturbed 
sinusoidal spin-spiral [Fig. f|(c)]. Most interesting is the 
case < H < H c , where the spin structure may be de- 
scribed as a soliton lattice: regions of the spin-flop phase 
are interrupted at regular intervals by magnetic domain 
walls, or solitons [Fig. [|(b)]. In each soliton the direction 
of staggered magnetization rotates by an angle tt. At 
H — > H c , the density of solitons starts to decrease very 
rapidly, as 1/ \n(H c — H)\ |J. The transition is thus al- 
most first order: as H — > H c , the two magnetic satellites 
at (1 ± C, ±C, 0) converge to eventually produce a single 
peak at (1, 0, 0). 

The exact expression for (3{H) [Eq. ([§] enables us to 
derive parametric equations for the field-dependence of 
magnetization and incommensurability parameter, and 
directly compare these predictions to experimental re- 
sults for Ba 2 CuGe 2 07. Using the formula M = —dF/dH 
and the equalities for derivatives of elliptic functions |T^] 
one gets: 
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The magnetization curve is continuous at the critical 
field, while x(H) = jjjj- diverges as H c is approached 
from below, and is constant and equal to x± at i? > H c . 
In zero field \ is equal to (x\\ + X_l)/2. We now use Eq. 
(^) together with the formula (||) to fit the experimen- 
tal x(H) for Ba 2 CuGe 2 07 measured at T = 2 K. With 
X± = 3.43 x 1(T 5 cmu/g, x\\ = 0.89 • 10" 5 emu/g and 
H c = 1.88 T a very good fit is obtained (Fig. g, solid 
line for -ff||c). In the refinement we have included a lin- 
ear contribution to x(-ff). This term is present in both 
Xc(H) and Xa(H) and is a result of non-linear corrections 
to the local susceptibility, neglected in our treatment. 

As H — > H c the soliton density decreases, i.e., the 
period of the magnetic structure increases with increasing 
field. In our case expressing £ in terms of /3 is rather 
straightforward and yields: 
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Again, Eq. (j^), when combined with Eq. (^), gives a 
parametric curve for Q{H), that can be fit to the ex- 



perimental data for Ba 2 CuGe 2 C>7 using H c as the only 
adjustable parameter [£(0) = 0.027 is measured indepen- 
dently]. The result with H c — 2.13 T is shown in a 
solid line in Fig. |[ A remarkable agreement is obtained 
except very close to H c . Discrepancies close to the tran- 
sition point are to be expected, since the transition is 
almost first order, the soliton lattice is very soft and eas- 
ily pinned by any impurities. The same effects prevent us 
from experimentally observing a true divergence in x(i?) 
at the critical field (Fig. ||). 

Finally, we can estimate the actual value of the critical 
field using the previously measured exchange constants 
and C(0). With J a b = 0.48 meV |ll| the exchange energy 
per bond is J = 2J a bS 2 « 0.24 meV. For x\\ and x± 
we can use the values for a classical Heisenberg antifer- 
romagnet at T = 0. ESR measurements (l^] provide us 
with the gyromagnetic ratios of Cu 2+ : g a = 2.044 and 
g c = 2.474. Using Eq. §) and X X = {9Vb) 2 /(8J) we 
obtain H c = 3.3 T, that should be compared to the ex- 
perimental value H c « 2.1 T. Considering that in these 
estimates we have completely ignored quantum- and tem- 
perature corrections to x and p Sl a 30% consistency is 
indeed acceptable. 

In summary, we have observed a rare type of CI tran- 
sition that is driven exclusively by the changing strength 
of the commensurate potential. The latter is directly 
controlled in an experiment by varying the magnetic 
field. A transition of this kind was envisioned over 
three decades ago by Dzyaloshinskii, and now we find 
that Ba 2 CuGe 2 C>7 exhibits it in its original form. Now 
that the underlying physics is rather well understood, 
Ba 2 CuGe 2 07 can be used as a very neat and simple 
model system for further studies of spiral magnetism. For 
example, it would be very interesting to look closer at the 
magnetic critical behavior. Another promising direction 
for future work is the study of the effect of an in-plane 
magnetic field. 
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FIG. 1. A schematic view of a Cu-Ge-O layer in 
Ba2CuGe2 07. The arrows indicate the components of 
Dzyaloshinskii vector D, allowed by the symmetry. 

FIG. 2. Field dependence of the magnetic susceptibility 
measured in Ba2CuGe2 07 at T = 2 K for the magnetic field 
applied along the c (circles) and a (crosses) axes of the crys- 
tal, respectively. The solid line is a theoretical fit to the data, 
as described in the text. 

FIG. 3. Field dependence of the magnetic propagation vec- 
tor in Ba2CuGe207 measured at T = 2.4 K. The solid line is a 
theoretical fit given by Eqs. (^Q). Insert: Elastic scans across 
the antiferromagnetic zone center measured in Ba2CuGe207 
at T = 2AK for two different values of magnetic field applied 
along the (0, 0, 1) direction. 

FIG. 4. Spin configurations for the spin-flop phase at 
H > H c (a), the soliton lattice at H = 0.997# c (b) and 
the circular spin spiral at H = (c). 
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